Plasmon and coupled plasmon-phonon modes in graphene 
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Plasmon and coupled plasmon-phonon modes in graphene are investigated theoretically within 
the diagrammatic self-consistent field theory. It shows that two plasmon modes and four coupled 
plasmon-phonon modes can be excited via intra- and inter-band transition channels. It is found 
that with increasing q and carrier density, the plasmon modes couple strongly with the optic- 
phonon modes in graphene. The coupled plasmon-phonon modes exhibit some interesting features 
which can be utilized to realize the plasmonic devices. Our results suggest that the carrier-phonon 
interaction should be considered to understand and explain the properties of elementary electronic 
excitations in graphene. 
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I. INTRODUCTION 

Since graphene was discovered in 2004 [l|, it has 
brought us many surprises and opened up a new field 
of research in condensed matter physics and nano- 
electronics. From a physics point of view, carriers in 
graphene are massless, gapless and relativistic Dirac 
particles with nearly linear energy dispersion. Thus 
graphene brings us a great opportunity to look into the 
many-body interactions in a Dirac quasi-particle sys- 
tem. Recently the elementary electronic excitation via 
plasmon modes resulting from Coulomb interaction has 
been intensively investigated for graphene based elec- 
tronic systems. Das Sarma and Hwang 0] and Wang 
and Chakraborty Q have found theoretically that the 
plasmon frequency of the Dirac fermions in graphene 
shows very different behaviors from the usual 2DEG sys- 
tems. More interestingly, it has been shown that the 
plasmon frequencies in graphene can be within the ter- 
ahertz (10 12 Hz or THz) bandwidth d and 1 - 10 THz 
sources based on plasmon amplification in graphene can 
be generated for moderate carrier densities by applying 
the gate voltages |5[. Moreover, Jablan et.al. has demon- 
strated that the plasmon excitation in graphene can have 
low losses for infrared frequencies @ and proposed that 
such a feature can be utilized for nano-photonic device 
applications. The results obtained from further theoret- 
ical investigations have indicated that the guided plas- 
mon waves in graphene p-n junctions can be achieved 
and graphene can then be applied as a material for nano- 
plasmonics Q- However, here we notice that most of the 
published theoretical studies mentioned above were car- 
ried out without considering the carrier-phonon coupling 
in graphene and only the intra-band transition channel 
was taken into consideration. Very recently, Bostwick 
et.al. observed the plasmarons in free-standing doped 



graphene, which results from bound states of charge car- 
riers with plasmons 0. 

On the other hand, recent experimental work has sug- 
gested that the adiabatic Born-Oppenheimer approxima- 
tion fails in graphene because of the presence of the in- 
teractions between the peculiar Dirac fermions and optic 
phonons [9]. It has been experimentally demonstrated 
that the plasmon modes in epitaxial graphene on SiC 
wafer can be strongly coupled with the surface optical 
phonon modes in silicon carbide [To| . Bostwick et.al. 
have found experimentally that the electron-electron in- 
teraction and the electron-phonon coupling should be 
considered on an equal footing to understand the proper- 



ties of Dirac quasi-particles in graphene [llj . Siegel et. al. 
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have experimentally demonstrated electron-electron in- 
teraction results in unique renormalizations, electron- 
phonon interaction diminish with decreasing doping in 
graphene [12j. The measurement of the Raman spec- 
trum has now been widely used to investigate the prop- 
erties of photo-excited electronic excitations in graphene 
systems [l3l [lij . The results obtained from recent Ra- 
man measurements have shown that the position of optic- 
phonon— induced G-peak depends very weakly on car- 
rier density in the sample. However, the Raman inten- 
sity around the G-peak increases with increasing doping 
concentration in graphene. Such an effect is an identi- 
fication that both the electron-electron interaction and 
the electron-phonon coupling affect strongly the elemen- 
tary electronic excitations in graphene. These experi- 
mental findings have clearly indicated that the carrier- 
phonon coupling is very important in understanding and 
explaining the electronic and optoelectronic properties 
of graphene, especially the elementary electronic excita- 
tions. 

We know that in conventional electronic gas systems 
the plasmon modes can be strongly modified by phonon 
scattering so that the new modes, namely the coupled 
plasmon-phonon modes, are formed. Raman scattering 
jl5j , infrared and photoemission measurements [l6j , and 
ultrafast pump-and-probe experiments [171 . etc. have 
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been widely applied to detect and study the coupled 
plasmon-phonon modes in semiconductor based bulk and 
low-dimensional electronic systems. The investigation 
of the plasmon and coupled plasmon-phonon excitations 
from electron gas systems is also a fundamental basis 
for the realization of these systems as plasmonic devices 
for various applications. Because graphene is an ideal 
2DEG system with high carrier density and high car- 
rier mobility, it is natural for us to expect that graphene 
can exhibit more excellent and unique plasmonic prop- 
erties for device applications. Thus, in order to gain an 
in-depth understanding of the many-body interactions in 
graphene systems and, on this basis, to explore the poten- 
tial applications of graphene as practical nano-plasmonic 
devices, it is necessary and significant to examine the 
roles which the carrier-carrier interaction and the carrier- 
phonon coupling can play in such a Dirac fermion system. 
And this becomes the prime motivation of the present 
theoretical investigation. 



II. THEORETICAL APPROACH 



the change of the carrier wavevector during a c-c scatter- 
ing event, V q = 2ire 2 / e^q is the two-dimensional Fourier 
transform of the Coulomb potentail, and 
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1 + aA 



kq 



5 a, 



(i-<W)> 



with A kq = (fc + <jcos6>)/|k + q|, B kq = (gsin0)/|k + q| 
and 9 being the angle between k and q. 

In the present study, we assume that the graphene sys- 
tem can be separated into the carriers of interest and the 
rest of graphene crystal. For the case of carrier interac- 
tions with 2D-like phonons, the interaction Hamiltonian 
takes a form 

H c _ p = W q a q e i(q ' r+a ^ + W*a q e-* (q ' r+a, ' i) , (3) 

where q = (q x , q y ) is the phonon wavevector along the xy- 
plane, (a q , a q ) are the canonical conjugate coordinates 
of the phonon system, Wq is the carrier-phonon (c-p) 
interaction coefficient, and uj q is the phonon frequency. 
The Fourier transform of the matrix element for bare c-p 
interaction can be written as 



Here we consider a graphene sheet in the sy-plane. 
A carrier (electron or hole) in a monolayer graphene 
can be described by the Dirac equation for a mass- 
less neutrino. The energy spectrum and wavefunction 
for a Dirac quasi-particle in graphene in the absence 
of the scattering mechanisms can be obtained analyti- 
cally. They are, respectively, E\(k) = Xy\k\ = Xjk and 
^ Ak (r) = |k,A >= 2- 1 / 2 [l,Ae^]e ik ' r in the form of a 
row matrix. Here, k is the wavevector for a carrier and 

k = ^Jk 2 + k 2 , r = [x, y), 7 = hvp is the band parameter 

with vf ~ 10 8 cm/s being the Fermi- velocity of a Dirac 
quasi-particle, A = +1 for an electron and A = — 1 for a 
hole, and <j> is the angle between k and the x direction. 

With the single-particle wavefunction for a carrier, we 
can calculate the electrostatic energy induced by bare 
carrier-carrier (c-c) interaction in graphene through 

VxfAiAJAa = S KXlKXz J d 2 T 1( l 2 r 2 ^ k /(ri)VA lkl (ri) 
xV(r 1 - T 2 )ipy 2 ^(r 2 )ipx 2 k 2 (r 2 ), 

(1) 

where V(r) = e 2 /e OC) |r| is the Coulomb potential with 
£00 ~ 1 being the high-frequency dielectric constant for 
graphene layer Q, and S Xl \ 2 \ 3 \ 4 = -i(Ai+A2+A 3 +A 4 )/2 = 
±1 or ±i is a sign function related to the charge sign of 
a carrier in different bands. After considering a momen- 
tum conservation law, namely the momentum that flows 
into the c-c scattering conserves with what flows out for 
different scattering processes, the bare c-c interaction be- 
comes 



V, 



a/3 



SapVqGap^k, q), 



(2) 



where we have defined a = (A'A) for intra- (i.e., A' = A) 
and inter-band (i.e., A' 7^ A) transition, q = (q x ,q y ) is 



A)H,tt)|t4(k,q)| 2 , 



(4) 



where f2 is the electronic excitation frequency, 
Do(uj q ,ft) — 2ftio q /[(ft£l) 2 — (hu q ) 2 ] is the bare phonon 
propagator, and \U\>\(\t, q)| 2 = | < k + q, A'| W q |k, A > 
| 2 . The published experimental results have indicated 
that the optic-phonon coupling is a major factor in 
determining the intensity and peak-position of phonon- 
related Raman spectrum in graphene fl3l [Til ]. Thus, 
in conjunction with these experimental findings, in 
this study we consider carrier interaction with only 
optic-phonons. On the basis of a valence-force-field 
model, the coupling coefficient for carrier interactions 
with long-wavelength optic-phonons in graphene is [l8j 



yv q 



where g = h^{B /b 2 ) / ^/2phwo with p ~ 6.5 x I0~ 8 g/cm 2 
being the areal density of the graphene sheet, lu = 196 
meV the optic-phonon frequency at the T-point, B ~ 
2 is a dimensionless parameter, and b = a/y/3 is the 
equilibrium bond length. Furthermore, 



M l = 







-i<j>„ 







and M q 







ie' 







, (6) 



for coupling with, respectively, the longitudinal (I) and 
transverse (t) phonon modes, where <p q is the angle be- 
tween q and the x-axis. Thus, the squares of the carrier- 
phonon scattering matrix elements are 



A ,(q, k)| 2 = («? 2 /2)[l - A'Acos(0 + 4>' - 24> q )\ 



and 



|t/* A ,(q,k)| 2 = {g 2 /2)[l + A'Acos^ + 0' - 20,)]. 
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The total contribution of the optic-phonon scattering is 

|f/ A v(q,k)| 2 = |C/ AA ,(q,k)| 2 + |C/AA<(q,k)| 2 = 5 V, (7) 

which is independent on q and k. This feature is very dis- 
tinct from that in the conventional semiconductor-based 
2DEG systems. 

In a diagrammatic self-consistent field theory [Tgi] , the 
effective c-c interaction in the presence of the c-p coupling 
is given by 

reff 



v:y = [v%+w*>]e-i, 



where 



toft = S a ^S(k) - [V£ c p -r ry 

dynamical dielectric function matrix element and 



W^ln^k^O) is the 



n A / A (k,q;fi) = g s g v 



fy[Ey(k + q)}~f x [E x (k)} 
Ml + E\, (k + q) - E x (k) + iS 



is the pair bubble or density-density correlation function 
in the absence of the c-c screening, with f\(x) being the 
Fermi-Dirac function for a carrier in the A band and g s = 
2 and g v — 2 counting respectively the spin and valley 
degeneracy. After summing e a p over k and setting j = 
(A'A) = 1 = (++), 2 - (+-), 3 = (-+) and 4 = (— ), 
the dielectric function matrix is obtained as 





e(g,0) 
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-a 2 1 - 


- a 3 ~b 3 
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- a 4 + b 4 _ 



(8) 



where 



a\>\ 



(q, fi) = -(V q /2) ^(1 + A'AA kq )n j (k, q; Q) 



is induced by the c-c interaction, and 

v (g, O) = -£>oK, A) V |C/ A 'A(k, q)| 2 n, (k, q; fi) 



?A'Al 



which results from intra- (i.e., j — 1 or 4) and inter-band 
(i.e., j = 2 or 3) electronic transitions. The modes of 
the elementary electronic excitation are determined by 
Re|e| -> 0. 

Now we consider a n-type (or positively gated) 
graphene with an electron density n e and a distribution 



function f e (x) = [1 + e 



(x-E F )/k 



bT ] 1 where Ep is the 



Fermi energy (or chemical potential) for electrons in the 
conduction band. In such a case, the valence band is fully 
occupied so that the distribution fh(x) = 0, the hole den- 
sity rih — 0, and 04 = 64 = 0. The approach can also 
be applied to the p-type (or negatively gated) graphene 
samples because the electrons and holes are symmetri- 
cal in graphene system. At a long-wavelength (q — > 0) 
and low-temperature (T — > 0) limit, the frequency of 
the collective excitation of the uncoupled plasmon mode 
(bj = 0) is obtained as 



J p = \j2e 2 E F q/e c 



q 



1/2 



,1/4 



(10) 



which is induced by intra-band transition within the con- 
duction band. Here, Ep = 7kp with k^= ^/7m e being 
the Fermi wavevector in graphene and q is defined as 
the plasmon wavevector or momentum. By solving the 
equation 



n 



AEp 



In 



2E F + n 



2E F - fl 



0. 



(11) 



is caused by the c-p coupling. In this study we use a 
matrix to present the dielectric function. In contrast to 
previous theoretical study in which only the intra-band 
transition is considered for calculating the dynamical di- 
electric function, the present work includes the contribu- 
tions from inter-band transition channels as well. Thus, 
the dielectric function matrix is a 4 x 4 matrix. In the 
absence of the phonon scattering, i.e., when bj = 0, the induced by intra-band transition for electrons within the 
dynamical dielectric function matrix becomes that ob- 
tained under the standard random phase approximation 
(RPA). The determinant of the dielectric function matrix 
Eq. (p} is given by 



we can obtain the uncoupled plasmon mode uj' p induced 
by inter-band transition for electrons from valence band 
to conduction band in graphene. Similar to a conven- 
tional 2DEG system, the single-particle (electron-hole 
pair) excitations (SPE) for intra- and inter-band transi- 
tion channels are also allowed in graphene. The intra- 
(upper) and inter- (lower) band transition boundaries 
are still given respectively by w s = 7|k + q| =F 7|k|. 
Thus, we are able to identify the intra-band SPE regime 
< lj s < jq for q < 2kp, and jq — 2Ep < lo s < jq 
for q > 2kp. The inter-band SPE regime is allowed for 
2Ep — "/q < lo s < 2Ep + 79 in graphene. Furthermore, 
a critical parameter q c =k^(2 + c — y/Ac + c 2 ) can be in- 
troduced to indicate the minimum momentum transfer 
with c = e 2 /(eoo7), where the plasmon mode enters to 
the electron-hole continuum. 

Including the electron-phonon coupling (bj ^ 0), we 
obtain two coupled plasmon-phonon modes, 

n ± = -i=K 2 + C 2 ± V /( Wo 2 - W 2)2 +WoW 3 ]l/2 (12) 



conduction band in graphene, where uj^ = 8g j q Ep/n 
Moreover, by solving the equation 

7rw 2 (fJ 2 - ujD + Ag 2 EpLu ^ 2 



= [(1 + ai + 6i)(l + a 4 + 64) - 0104] 
x [(1 - a 2 - & 2 )(1 - 03 - 63) - a 2 a 3 ] 



2E f 



(9) 



~ 4irE F {n 2 -uj 2 ) + 16g 2 E 2 F ujo 2Ep — Q ~ ' 

we acquire two coupled plasmon-phonon modes and 
f2 2 induced via inter-band transition for electrons from 
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FIG. 1: Dispersion relation for uncoupled plasmon modes 
uj p (solid line) and lj' p (dashed line) at a fixed electron den- 
sity n e = 1 x 10 12 cm -2 . The regime III (lightly shaded) 
and regime IV (heavily shaded) correspond respectively to 
the inter-band SPE and intra-band SPE. The regime I and 
regime II are for undamped plasmon excitation and the min- 
imum momentum transfer q c ~ 0.51ki? is indicated with 
k F ~ 1.77 x 10 6 cm' 1 . 



valence band to conduction band in graphene system. 
The analytical results presented in this Section show that 
after including the inter-band transition channels and the 
electron-phonon coupling, the new modes of elementary 
electronic excitations can be achieved in graphene. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 



In Fig. 1, the dispersion relation for uncoupled plas- 
mon modes, uj p and cu' p induced respectively by intra- 
and inter-band transitions, are shown at a fixed electron 
density in graphene. As being pointed out by other au- 
thors 2, 3], the plasmon frequency induced by intra-band 

excitation in graphene, lo p ~ q 1 ! 2 ~ nl^ 4 , is acoustic- 
like and depends strongly on q. The density dependence 
w p ~ n~J 4 in graphene differs from that lu p ~ nj 2 in a 
conventional 2DEG system. It shows that the plasmon 
frequency u p of graphene can be effectively controlled 
through tuning carrier densities by gated voltages. In the 
presence of the inter-band transition channels, we obtain 
a new uncoupled plasmon mode uj' p induced by excitation 
of electrons from occupied valance band to the empty 
states in the conduction band. Due to Pauli exclusion 
principle, uj' p > 2Ep is optic-like and depends relatively 
weakly on q. With increasing q, uj p and lu' p converge in 
regime III. We find that uj p depends more strongly on car- 
rier density and q than to' p does. For a undoped graphene 
sample so that n e —¥ 0, the plasmon frequency uj p —¥ 
and the excitation is restricted, whereas there exits the 
plasmon mode with a frequency w' due to inter-band 
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FIG. 2: Dispersion relation for coupled plasmon-phonon 
modes fi+ (upper panel) and f2_ (lower panel) induced 
by intra-band transition for different electron densities in 
graphene as indicated. Here uio — 196 meV is the optic- 
phonon frequency. 
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FIG. 3: Dispersion relation for coupled plasmon-phonon 
modes fii (~ uuo) and Q.2 (> 2Ef) induced by inter-band 
transitions for different electron densities n e — 1 x 10 12 cm~ 2 
(solid line), n e = 3X10 11 cm -2 (dashed line) and n e = lxlO 11 
cm -2 (dotted line). The three curves for Q.\ coincide roughly. 



electronic transition. Thus, this new plasmon mode can 
be used to investigate the collective electronic excitation 
from an nearly intrinsic graphene. In Fig. 1, we also show 
the regimes III and IV for q <kp within which the SPE 
can take place. Similar to a conventional 2DEG, within 
the SPE regime the Landau damping occurs so that the 
imaginary part of the dielectric function Ime(u, q) =/= 0, 
and that if the plasmon excitation exists in these regimes 
for q > q c ~ 0.51k F (see e.g. Ill and IV in Fig. 1), 
the plasmon modes are strongly damped and decay into 
the electron-hole pairs. Thus, one can see that the plas- 
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mon mode w p is undamped for q < q c in regime I where 
Ime(u;, q) — 0, whereas the plasmon mode u' p is damped 
in the full regime of III in which Ime(u),q) 7^ 0. Conse- 
quently, in graphene, the undamped plasmon excitation 
can be achieved via intra-band transition, while the plas- 
mon excitation via inter-band transition is damped due 
to the Landau damping effect. 

In Fig. 2, the dispersion relation for coupled plasmon- 
phonon modes f2+ and induced by intra-band tran- 
sition is shown for different electron densities. For small 
values of q < qo, the lower frequency branch f2_ < luq is 
acoustic-like (or plasmon-like) which depends strongly on 
q, whereas the higher frequency branch f2+ ~ ujq is optic- 
like (or phonon-like) which depends rather weakly on q. 
Interestingly, for relatively large values of q > q a , the 
lower frequency branch becomes optic-like and f2_ ~ u)q, 
whereas the higher frequency branch is acoustic-like and 
f2+ > ujq. When q > qo, the plasmon modes couples 
strongly to the optic-phonon modes and the frequencies 
of the coupled plasmon-phonon modes in graphene differ 
significantly from the corresponding plasmon and optic- 
phonon frequencies. The qo value, at where the two 
branches of the excitations change the nature of the q 
dependence, decreases with increasing the electron den- 
sity. It is interesting to compare the results for the cou- 
pled plasmon-phonon modes in graphene with those in a 
conventional 2DEG (C2DEG) system. In a semiconduc- 
tor based C2DEG, there exit three modes for coupled 
plasmon-phonon excitations induced by intra-subband 
transition channels via electron interaction with optic- 
phonons through Frohlich coupling [20]. They are: i) 
longitudinal-optic (LO) phonon-like mode w+ ~ ^lo 
with lolo being the LO-phonon frequency; ii) transverse- 
optic (TO) phonon-like mode oj_ ~ luto with wto being 
the TO-phonon frequency; and hi) plasmon-like mode 
uji, respectively. In a C2DEG, the plasmon-like mode 
does not coupled markedly with the phonon-like modes 
and the phonon-like modes depend quite weakly on q 
[2fj| . The reason why only two coupled plasmon-phonon 
modes are observed in graphene is that at the T-point, 
ujlo — uto — w o f° r graphene [l8[ in contrast to a 
C2DEG where normally ui^o 7^ ^>to at the T-point. The 
different q-dependence of the coupled plasmon-phonon 
modes in graphene and in a C2DEG is mainly induced 
by different energy spectra in two electronic systems. As 
we know, the elementary electronic excitation is achieved 
via varying the energy and momentum of the electrons. 
Thus, the energy spectrum of the electronic system de- 
termines mainly the dispersion relation of the electronic 
excitation. In contrast to a parabolic energy spectrum 
for a C2DEG, graphene has a linear energy spectrum. As 
a result, the dispersion relation of the coupled plasmon- 
phonon modes in graphene differs significantly from those 
in a C2DEG. From Fig. 2, we can also hnd that with 
increasing electron density n e , the go value decreases, 
and the plasmon and phonon modes are strongly cou- 
pled. This is because the plasmon frequency uj p ~ n~J^ 
in graphene can be comparable to the optic-phonon fre- 



quency ujq at a relatively large n e . 

In Fig. 3, we show the dispersion relation for cou- 
pled plasmon-phonon modes with frequencies fii and Q2 
induced by inter-band transition in graphene for differ- 
ent electron densities n e . We see that a strong optic- 
phonon— like mode with a frequency f2i ~ ujq can be 
generated in graphene and the frequency of this mode 
depends very weakly on both q and the electron density. 
In contrast, the plasmon-like mode ^2 ~ u' p depends sen- 
sitively on electron density and on plasmon wavevector 
q. 0,2 increases with n e due to its plasmon-like nature of 
the excitation. We note that both modes with f2i and 
Q2 are optic-like and uncoupled with each other. For 
a n-type graphene, an excitation energy E > Ep is re- 
quired to excite electrons in the valence band into the 
conduction band and to induce the elementary electronic 
excitations due to the phase-space restriction. As a re- 
sult, there are always f2i ~ ujq and VL2 > 2Ep for coupled 
plasmon-phonon frequencies induced by inter-band tran- 
sition in graphene and these modes are undamped by the 
Landau-damping effect. 

It is known that the elementary electronic excitation 
can be achieved in an electron gas system through elec- 
tronic transition from occupied lower-energy states to the 
empty higher-energy states. This includes both intra- 
and inter-band excitations. In a semiconductor based 
C2DEG, because the band-gap between the conduction 
band and the valence band is quite large, the inter-band 
electronic transition requires a quite large excitation en- 
ergy and thus is less possible. However, graphene is a 
gapless 2DEG and the inter-band electronic transition in 
graphene can therefore be achieved via a small energy 
transfer. Hence, the elementary electronic excitation via 
plasmon and coupled plasmon-phonon modes induced by 
inter-band transition can be more possibly achieved in 
graphene than in a C2DEG. Because the carrier density 
in graphene can be effectively modulated up to ~ 10 13 
cm~ 2 by applying the gate voltages H, H3, EH, the plas- 
mon and coupled plasmon-phonon excitations can be ex- 
perimentally detected and studied by different methods 
[ToL [l7| . It is found experimentally that the position of 
the G-peak induced by optic-phonon scattering in the Ra- 
man spectrum in graphene varies very little with varying 
the carrier density in the sample. However, the Raman 
intensity around the G-peak depends quite strongly on 
the carrier density in graphene. The results from the 
present study show that for small values of q, optic- 
phonon— like modes can be generated through coupled 
plasmon-phonon excitations via intra-band (see f2_|_ in 
Fig. 2 for q < qo) and inter-band (see Q.% in Fig. 3) elec- 
tronic transition channels. Such modes depend rather 
weakly on q and n e for small values of q. This suggests 
that the position the optic-phonon— induced G-peak in 
the Raman spectrum does not change significantly with 
varying the carrier density in the sample, in line with 
the experimental finding. From Fig. 2, we see that for 
relatively large values of q > qo which decreases with 
increasing n e , the frequency of optic-phonon— like mode 
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induced by intra-band excitation depends strongly 
on electron density. This implies that the Raman in- 
tensity around the G-peak can be altered with vary- 
ing the carrier density in graphene. Because the optic- 
phonon— like mode induced by inter-band excitation fl\ 
depends very little on n e (see Fig. 3), the change of the 
Raman intensity around the G-peak is mainly induced 
by intra-band excitation mechanism. Thus, our results 
can be applied to understand and explain the experimen- 
tal findings. These experimental and theoretical results 
indicate that the electron-phonon interaction can affect 
strongly the features of elementary electronic excitations 
in graphene. Furthermore, the results shown in Figs. 2 
and 3 demonstrate that the coupled plasmon-phonon ex- 
citations in graphene can provide frequency-tunable plas- 
monic modes. Hence, graphene can be used as plasmonic 
device for various applications. 

IV. CONCLUSIONS 

In this work, we have developed a tractable theory 
approach to study the plasmon and coupled plasmon- 
phonon modes in graphene in which the carrier-carrier in- 
teraction and the carrier-phonon coupling are taken into 
consideration. We have also included the effect of the 
presence of the inter-band electronic transition channels 
within our calculations. The main conclusions obtained 
from this study are summarized as follows. 

An acoustic-like plasmon mode with a frequency uj p ~ 

q 1 / 2 nl^ 4 can be generated via intra-band electronic 
transition in graphene and the carrier density dependence 
of this mode differs from that in a conventional 2DEG. 
This is in line with previous theoretical finding. We have 
found that an optical-like plasmon mode with a frequency 



uj' p can be excited via inter-band electronic transition in 
graphene. However, this mode is within the Landau- 
damping regime. 

We have found that four coupled plasmon-phonon 
modes can be generated in graphene. The features of 
two coupled-phonon modes induced by intra-band tran- 
sition channels differ significantly from those observed 
in a conventional 2DEG system. Two optic-like modes 
for coupled plasmon-phonon excitation can be observed 
via inter-band electronic transition channels. Three cou- 
pled plasmon-phonon modes depends strongly on elec- 
tron density and plasmon wavevector. 

The results obtained from this study demonstrate that 
if the inter-band transition channels are taken into con- 
sideration, new modes for plasmon and coupled plasmon- 
phonon excitations can be obtained in graphene. The 
properties of the plasmon and coupled plasmon-phonon 
excitations in graphene are very distinct from those 
in the conventional 2DEG systems. Our results con- 
firm that the carrier-carrier interaction and the carrier- 
phonon coupling should be equally considered to under- 
stand the properties of elementary electronic excitation 
in graphene. Finally, our results presented and discussed 
in this paper can be used to understand and explain the 
experimental findings obtained from, e.g., the Raman 
measurements. 
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